Surface-induced charge state conversion of nitrogen-vacancy defects in nanodiamonds 
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We present a study of the charge state conversion of single nitrogen- vacancy (NV) defects hosted 
in nanodiamonds (NDs). We first show that the proportion of negatively- charged NV - defects, with 
respect to its neutral counterpart NV°, decreases with the size of the ND. We then propose a simple 
model based on a layer of electron traps located at the ND surface which is in good agreement with 
the recorded statistics. By using thermal oxidation to remove the shell of amorphous carbon around 
the NDs, we demonstrate a significant increase of the proportion of NV - defects in 10-nm NDs. 
These results are invaluable for further understanding, control and use of the unique properties of 
negatively-charged NV defects in diamond. 

PACS numbers: 76.67.-n, 81.05.uj, 81.07.-b, 71.55.-i 



The negatively-charged nitrogen- vacancy (NV - ) 
defect in diamond has received considerable interest 
over the last decade owing to its perfect photostability 
and to its ground-state electron spin properties, which 
combine a long coherence time^ and the ability to 
undergo spin- sensitive optical transitions under ambient 
conditions^. Such properties are at the heart of a wide 
range of emerging technologies, from imaging in life 

urn 



7 to high-resolution magnetic sensing 7 y and 

quantum information processings^. For most of these 
applications, it is essential to use NV - defects as close 
as possible to the diamond surface, either to enhance 
the sensitivity of diamond-based magnetic sensing or 
to efficiently couple the NV - defect photoluminescence 
to a photonic waveguide or a microcavit ^ 12 * 13 *. This 
condition can be fulfilled by co nsiderin g NV - centers 
hosted in a nanodiamond (ND)^^^ or created by 
shallow implantation in a single-crystal bulk diamond 16 . 
However, it was recently reported that the preferred 
charged state of NV defects close to the diamond 
surface is the neutral stat d 16 * 1 ^ NV°, which does not 
exhibit the appealing optical and spin properties of 
its negatively-charged counterpart. It is then highly 
desirable to understand the processes leading to charge 
state conversion, with the aim of producing efficiently 
NV - color centers close to a diamond surface. 

In this context, we investigate the charge state 
conversion of single NV defects hosted in NDs of variable 
size. Since NV° and NV - defects can be discriminated 
through their photoluminescence spectrum (Fig. TJ, we 
determine the proportion V^y- of negatively-charged 
NV - defects found in an ensemble of NDs containing 
individual NV color centers, as a function of the ND size. 
The measurements are done with a confocal microscope 
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Figure 1: (color online) Photoluminescence (PL) spectra (nor- 
malized to their respective maximum value) of single NV - 
(red curve) and NV° (blue curve) color centers in nanodi- 
amonds. The zero-phonon line (★ symbols) of NV - (resp. 
NV°) emission is located at 637 nm (resp. 575 nm). The in- 
set shows the atomic structure of the NV defect, consisting of 
a substitional nitrogen atom (N) associated to a vacancy (V) 
in an adjacent lattice site of the diamond crystalline matrix. 



combined with an atomic force microscope. We show 
that V NV - decreases with the ND size, indicating that 
surface-related effects play a significant role in charge 
state conversion between NV - and NV°. Many studies 
have indeed established that various impurities at a dia- 
mond surface can strongly modify the charge transport 
properties near the surface 18 . In the specific case of irra- 
diated and annealed diamond, graphitic related defects 
formed on the diamond surface are possible electron trap 
candidates 19 which can ionize donors required for an 
efficient conversion of the NV defects into the negatively- 
charged stated Following RefP^, a simple scheme is 
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then proposed for modeling the experimental observa- 
tions. By removing graphitic related defects through 
thermal oxidation^!, we finally show that the proportion 
of NV - defects relative to their neutral NV° charge state 
can be strongly enhanced, especially in the smallest NDs. 

We started from commercially available diamond 
nanocrystals with different particle size distributions 
(SYP 0.05 and 0.1, Van Moppes SA, Geneva). Such 
nanocrystals are produced by milling type-lb high- 
pressure high-temperature diamond crystals with a high 
nitrogen content ([N] « 200 ppm). The formation of 
NV defects was carried out using high energy (13.6 MeV) 
electron irradiation 21 followed by annealing at 800° C un- 
der vacuum during two hours. The electron irradiation 
creates vacancies in the diamond matrix while the an- 
nealing procedure activates the migration of vacancies 
to intrinsic nitrogen impurities, leading to NV defect 
bonding 22 !^. After annealing, the irradiated NDs were 
cleaned using a "piranha" solution which is a strong oxi- 
dizer intended to purify the ND surface 24 . After sonica- 
tion and washing with distilled water, aqueous colloidal 
solutions of dispersed NDs were obtained. A few drops of 
ammonium hydroxide were then added to the solutions 
to keep them stable for several weeks. NDs were finally 
deposited by spin-coating onto a piranha-cleaned silica 
coverslip. 

With the aim of studying the effect of the ND size on 
the charged state of NV defects, two solutions of irradi- 
ated NDs with different size distributions were prepared. 
These solutions correspond to a respective mean size of 
35 nm (SYP 0.05) and 65 nm (SYP 0.1), as determined 
by dynamic light scattering measurements. Due to the 
broad size distribution, NV defects in NDs with size rang- 
ing from 10 to 100 nm could be investigated with these 
two primary solutions. 

The photoluminescence properties of NV defects 
hosted in NDs were studied using a scanning confocal 
microscope integrated with an atomic force microscope 
( AFM) , as shown in Fig. [2] A laser operating at 532 nm 
wavelength was tightly focused onto the sample through 
a high numerical aperture microscope objective (x60, 
NA = 1.35). The photoluminescence emitted by NV 
defects was collected by the same objective and spec- 
trally filtered from the remaining pump light. Following 
standard confocal detection scheme, the collected light 
was then focused onto a 50 /xm diameter pinhole and di- 
rected either to a photon counting detection system for 
light intensity measurement or to an imaging spectro- 
graph allowing to distinguish between the NV - and NV° 
emission (see Fig. [I]). In addition, the AFM was used to 
record topography images of the sample, thus providing a 
direct in situ measurement of the ND size. For each pho- 
toluminescent ND, the emission spectrum and the ND 
size were measured, allowing to correlate the size to the 
charge state of the NV defects. 

In order to avoid possible bias in the statistics, the 
whole study was restricted to NDs hosting single NV 
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Figure 2: (color online) (a)-Sketch of the experimental setup 
combining an atomic force microscope (AFM, Asylum Re- 
search MFP-3D-BIO) with an inverted confocal microscope. 
NDs: nanodiamonds; O: oil immersion microscope objective 
(Olympus, x60, NA = 1.35); DM: dichroic beamsplitter; F: 
580 nm long-pass filter; PH: 50 fim diameter pinhole; FM: 
flip mirror directing the collected photoluminescence either 
to an imaging spectrograph equipped with a back-illuminated 
cooled CCD array, or to a Hanbury Brown and Twiss interfer- 
ometer consisting of two silicon avalanche photodiodes (APD) 
placed on the output ports of a 50/50 beamsplitter. (b),(c)- 
Photoluminescence raster scan of the sample (c) with the cor- 
responding topography image (b) . Red dotted circles indicate 
NDs hosting NV defects. The correspondence between the 
two images can be easily inferred. (d)-Typical second-order 
autocorrelation function g^ 2 \r) recorded for a single NV de- 
fect in a diamond nanocrystal. 



defect. The unicity of the emitter was systematically 
checked using photon correlation measurements. Since 
after the emission of a first photon, it takes a finite time 
for a single emitter to be excited again and then emit a 
second photon, a dip around zero delay r = appears in 
the normalized second-order correlation function 
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where I(t) is the photoluminescence intensity at time t. 
The later function is deduced from the histogram of time 
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Figure 3: (color online) Proportion of negatively- charged 
NV - defects V^y- as a function of the size of the photo- 
luminescent ND. The statistics is determined from a set of 
146 single NV defects and the size binning in the histogram 
was adapted in order to have an approximately equal number 
of investigated color centers per size slot. The red solid line 
is a fit of the data using equation 



delays r between two consecutively detected photons, 
recorded by using two avalanche photodiodes placed on 
the output ports of a beamsplitter, following the stan- 
dard Hanbury Brown and Twiss scheme (see Fig. [2|a)). 
The observation of an anticorrelation effect at zero time 
delay, g^ 2 \0) < 0.5, is the signature that a single NV 
defect is addressed (see Fig. [2|d))^. 

For a set of single NV defects, the emission spec- 
trum and the ND size were measured. Single NV° 
and NV - defects were then sorted according to their 
photoluminescence spectrum, in order to infer the 
proportion V^y- of negatively-charged NV - defect as a 
function of the ND size. The results of this experiment 
give a clear evidence that V^y- drastically decreases 
with the ND size (see Fig. [3|. Indeed, for 100-nm 
nanocrystals V^y- is close to unity but decreases to 
40% for 20 nm NDs. We noted that no blinking was 
observed in the photoluminescence intensity for all 
studied NV defects. In addition, all the measurements 
were performed while exciting NV defects belo w the 
saturation, in order to avoid photo-ionization effect <p6127| 



In the following, we introduce a simple model in order 
to explain the experimental observations. After the elec- 
tron irradiation, the vacancies migrate upon annealing at 
800° C and can be trapped by a substitutional nitrogen 
atom, thus forming a neutral NV° defect according to 
the reaction 



N + V -> NV° . 



(2) 



In order to convert the defect into its negatively- charged 
state NV - , it is widely accepted that an additional elec- 
tron donor located in the vicinity of the NV defect is 



requirecP^H. Since the most abundant electron donors 
in type lb diamond are substitutional nitrogen atoms, 
the process leading to the formation of NV - defect in- 
volves a second nitrogen atom and can be described by 
the "chemical" reaction^ 



NV° + N <=► NV~ 



• N 



(3) 



which can be characterized by an equilibrium constant 
K defined as 



K = 



[NV-] eg X [N+] e 
[NV°] e9 X [N] eg 



(4) 



where 



represents the volumic concentration of the 



chemical species i inside the ND. 

In addition, we consider that extra electrons provided 
by the substitutional nitrogen atoms can be captured by 
electron traps T located at the ND surface. This an- 
tagonist effect can be due to graphitic defects associated 
to a sp 2 -hybridized carbon layer created at the ND sur- 
face during electron irradiation and annealing 19 . Conse- 
quently, a part of nitrogen impurities can be ionized by 
this electronic depletion layer according to the reaction 



N + T -> N + 



(5) 



which prevents the nitrogen impurities to donate the re- 
quired electron allowing to convert the NV defect into its 
negatively- charged state, as represented by equation 

In the following, we assume that the reaction ([5| is to- 
tal. Furthermore, the content of NV defects per unit 
volume is supposed to be much smaller than all the 
other chemical species. This assumption is verified in 
our experimental conditions since the study was limited 
to NDs hosting single NV defects. When the initial vo- 
lumic concent of traps is such that [T]^ > [N], all the 
nitrogen atoms are ionized and the charge state conver- 
sion characterized by reaction pi becomes fully ineffi- 
cient. On the other hand if [T]i < [N], the content of 
ionized nitrogen at equilibrium [N + ] e(? is fixed by the re- 
action (p5|, i.e. [N+] eg = [T]i and [N] eq = C - [T] i? where 
Co = [Nje g + [N + ]eg is the total nitrogen concentration. 

Since the electron traps are assumed to be located at 
the ND surface, their equivalent volumic content can be 
written as [T]i = cr^/a, where a is the characteristic size 
of the ND and ax the trap surface density. Using equa- 
tion Q , the proportion of negatively-charged NV defects 
V^y- is then given by 



NV- 



[NV~ 



1 



1 



[nv-u + [Nv°] eg 1+K (j* i - 1 y 

(6) 

where a t = cr T /Co is a threshold size. If a ^> a t , 
the number of nitrogen atoms is much higher than 
the number of electron traps, leading to V^y- = 1- 
Conversely, in the limit where a = at, i.e. Co = [T]i, 
all the nitrogen atoms are ionized by the surface traps, 
thus preventing the formation of any negatively-charged 
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NV defect in the nanocrystal and V NV - = 0. Using 
equation ([6| as a fitting function, a good agreement with 
experimental data points is achieved for K = 1.0 ± 0.5 
and a t — 14 ± 1 nm (see Fig. [3|. Considering a 200-ppm 
nitrogen content, which corresponds approximatively 
to 0.04 nitrogen atoms per cubic nanometer, the result 
of the fit leads to a surface density of electron traps 
<jt ~ 0.5 nm -2 . 

As already mentioned, graphitic related defects formed 
on the nanodiamond surface during electron irradiation 
and annealing are possible electron trap candidates^. 
In order to check the presence of graphitic shells or 
amorphous carbon layer, the composition of the studied 
NDs was analysed by recording their Raman spectra. As 
shown in Fig. [4|a) (red line), the sp 2 -hybridized carbon 
is evidenced through a broad D-line at 1350 cm -1 and 
a broad G-band around 1600 cm -1 , while the diamond 
core gives the usual sharp Raman line at 1335 cm~®. 
This result indicates that a significant content of sp 2 
carbon was partly covering the ND surface, even after 
acid cleaning with the piranha solution. 

According to the previously developed model, the 
proportion of negatively-charged NV defect in NDs 
might be improved by reducing the surface density <tt 
of electron traps. For that purpose, the graphitic shell 
which could be associated to the electronic depletion 
layer at the ND surface can be efficiently removed 
through surface oxidation. This can be achieved with a 
repetitive cleaning with a strong acid treatment, using 
a mixture of sulfuric, nitric and perchloric acid d 4 * 14 *. 
Another method, which does not require the use of 
toxic or a ggressive chemicals, consists in heating NDs 
in ai j 4 | 32 | 33 i This procedure leads to oxidation of 
sp 2 - hybridized carbon (graphite and disordered carbon), 
thus resulting in a higher purity of sp 3 -bonded carbon 
(diamond) 34 . In the following, we check if the propor- 
tion of negatively-charged NV - defect is modified after 
elimination of sp 2 -bonded carbon with air oxidation. 
For that purpose, the same colloidal solutions of NDs 
were oxidized in air at 550° C during two hours. We 
note that the stability of the ND aqueous suspensions 
remained unchanged after this surface oxidation. 
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Figure 4: (color online) (a) Raman spectra associated with 
NDs cleaned with a piranha solution after electron irradia- 
tion and annealing (upper red curve) and with subsequent air 
oxidation at 550° C during one hour (bottom blue curve). The 
G-band at 1600 cm -1 and the D-band between 1250 cm -1 and 
1450 cm -1 are associated to sp 2 -hybridized carbon atoms. Air 
oxidation leads to a higher purity of sp 3 -bonded carbon since 
both the D-band and the G-band are suppressed. Insets show 
the pictures of the aqueous colloidal solutions before and af- 
ter air oxidation. (b)-Proportion of negatively- charged NV - 
defects V NV - as a function of the ND size after air oxida- 
tion (blue points). The statistics is built from a set of 183 
single NV centers. The blue solid line is a fit obtained using 
equation The red points are reproduced from Fig [3] for 
comparison. 



The efficiency of the surface oxidation was first 
analysed by recording the Raman spectrum of the 
oxidized NDs. As shown in Fig. 4^a) (blue curve), 
the spectral features related to sp 2 -bonded carbon are 
greatly suppressed. The effect of surface oxidation was 
also confirmed by the color of the solution of NDs, 
which became whitish after air oxidation (see insets 
in Fig. [I|a)) 32 . In addition, the surface oxidation led 
to a size reduction of the NDs 35 . Using dynamic light 
scattering, the mean value of the size distribution was 
found to be reduced from 35 nm (resp. 65 nm) to 27 nm 
(resp. 53 nm) for the solution of SYP 0.05 (resp. SYP 
0.1). 

From a set of single NV defects in oxidized NDs, 



we inferred the proportion V^y- °f negatively-charged 
NV - defects as a function of the ND size (Fig. Qb)). 
As previously observed, V N y- decreases with the ND 
size. However, the proportion of NV - defects found 
in small NDs is strongly enhanced compared to the 
non-oxidized ND solution. Indeed, V^y- reaches 70% 
for near 10-nm NDs, which is two times larger than the 
result obtained for non-oxidized NDs (see Fig. |3|. By 
fitting the data with equation (J6|, a good agreement is 
found with parameters K = 1.4 ± 0.4 and a t = 6 ± 1 nm 
which corresponds to a surface density of charge traps 
<tt ~ 0.25 nm -2 . This value is twice smaller than 
the one obtained for non-oxidized NDs. According to 
modifications in the Raman spectra, we assume the 
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existence of other trapping sites inside the nanocrystal 
volume, being for instance associated to dislocations of 
the diamond matrix or to extended defect s formed by 
vacancies during the annealing process^^H 

Finally, we note that the photoluminescence intensity 
of single NV defects in oxidized NDs was improved by 
roughly 4-fold and the excited-state lifetime was found 
longer, as recently reported in the context of ensemble 
measurements^. Both effects can be understood by 
considering the layer of sp 2 - hybridized carbon atoms as 
a photoluminescence quencher. However, a quantitative 
study of the modifications of the photoluminescence 
properties of the NDs can be hardly done since both the 
emission intensity and the excited-state lifetime depend 
on the dipole orientation of the NV defect which is 
randomly distributed over each ND. 

To summarize, we have reported a study of surface- 
induced charge state conversion of single NV defects 
hosted in NDs. This study illustrates the influence of 
surface-related defects associated with sp 2 -hybridized 
carbon on the charge state of the NV color center. A 



simple model led to a threshold size for NDs hosting 
negatively- charged NV - defects, which is directly related 
to the surface density of electron traps. By using a ther- 
mal oxidation allowing to remove the shell of amorphous 
carbon around the NDs, we demonstrated that this 
threshold size can be significantly reduced, thus leading 
to an enhancement of the proportion of negatively- 
charged NV - defect in NDs of near 10 nm size. These 
results are invaluable for further understanding, control 
and use of the unique properties of negatively-charged 
NV defects in diamond, from high-resolution magnetic 
sensing to bioimaging. 
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